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Abstract 

We derive a method for a common treatment of Z' exchange, QED corrections, and weak 
loops. It is based on the form factor approach to the description of weak loop corrections 
to partial Z widths and cross sections. Problems connected with Z,Z' mixing are 
discussed with special care. Our theoretical results are applied to the package ZFITTER. 
We demonstrate two different ways to the data analysis - one based on an extension of 
the standard model cross sections, the other on model-independent formulae together 
with the Z width calculations in presence of a Z'. 

With the resulting package ZEFIT0ZFITTER, LEPl data from fermion pair production, 
including Bhabha scattering, can be analysed on, but also off the Z peak. Further, the 
code may be used at very high energies, e.g. in the region of a possible future linear 
e'''e“-collider. 
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1 Introduction 


Although the Standard Model [|^ has been verified with a precision including one loop cor¬ 
rections, there is a general consensus that we are far away from a final understanding of the 
elementary particle world. A unification of forces seems to happen at much higher mass scales 
than are directly accessible by present accelerators. Candidates for a truly unifying theory 
usually predict additional gauge degrees of freedom, thus leading in a natural way to the exis¬ 
tence of new, heavy neutral gauge bosons besides the photon and the Z boson of the standard 
theory (see e.g. p). 

Since 7 , Z, Z' are neutral particles with vector and axial-vector couplings, a search for a Z' is 
complicated by the fact that there is no special final state signature. Below the production 
threshold, theoretical predictions consist of minor quantitative modifications of the neutral 
current cross sections. As a consequence, one needs very precise predictions for cross sections 
and asymmetries. Important reactions for a dedicated search are: 

e+e-^(7,Z,Z0 —/+/-(7), (1) 

e+e“ —> ( 7 , Z, Z') —> e+e“( 7 ), (2) 

ep-^( 7 ,Z,Z')^eX( 7 ). (3) 

In earlier studies, we investigated QED (i.e. pure photonic) corrections together with Z' 
exchange both in e^e“-annihiIation into fermions [Q and in ep-scattering [^. In the latter 
reaction, we applied also part of the material presented here iZ' effects in presence of weak 
loops without Z,Z' mixing). 

In principle, Z' effects can be searched for by three different effects: 

• via virtual Z' exchange (at sufficient energy, also present without Z,Z' mixing) 

• via modification of the mass of the standard Z boson seen at LEPl due to Z, Z' mixing 

• via modification of the couplings of the standard Z boson seen at LEPl due to Z, Z' mixing; 
this in fact concerns two different, although related observables - the Z width [~ peak height] 
and cross sections [~ line shape]. 

For large enough Z' masses, the direct cross section contributions due to Z' exchange may be 
neglected at LEPl energies completely. Then, one can concentrate on the consequences of the 
Z,Z' mixing. In fact, LEPl is the ideal place to seareh for this phenomenon. 

From existing measurements it is known that the mixing is very small if not vanishing; see e.g. 
[Q and references quoted therein. In such a situation, one has to study carefully the interplay 
of the weak standard theory loop effects with the Z' infiuence via particle mixing. 

Often the LEPl data are analysed after a model-independent interpretation of the line shape 
(and of asymmetries) in terms of e.g. partial widths (or effective couplings), confronting those 
as ’observed quantities’ with theoretical expectations. Our approach will also allow to do so. 
In addition, we derived formulae which open the possibility of a direct interpretation of e.g. 
the line shape as function of the energy in terms of standard model parameters and, in parallel, 
of the Z, Z' mixing. (Not mentioning for a moment the direct Z' cross section part.) This 
is done in a specific way following the form factor notation of weak loop effects on which the 
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Dubna/Zeuthen approach (and e.g. the code ZFITTER which has been supplemented here) so 
heavily relies. 


We remind the reader that for each scattering process, there are four complex-valued form 
factors p, Ke, ^/, Kef, which allow an exact description of the weak (non-photonic) contributions 
to the scattering process M- For an exact description of Z,Z' mixing, one has to understand 
how these form factors become modihed. 

Since the partial and total widths of the Z boson [§] contain similar form factors pz,Kz, a 
similar procedure has to be applied there, too. 

We will also give a short comment on the relation of the formulae derived here to those used 
by other authors. 

This concerns the change of form factors by the Z,Z' mixing. 

Further, we dehne the weak mixing angle sin^ 6w from the unmixed mass parameters also in 
presence of Z,Z' mixing. As a result, the calculation of the standard weak loop corrections 
remains unchanged, too. Many other authors prefer to redehne the weak mixing angle as being 
derived from the physical mass of the Z mass eigenstate [^. This leads to the introduction 
of a p parameter in the dehnition of the weak mixing angle. The difference between the two 
procedures is essentially a different book keeping and does not influence the experimental 
determination of e.g. the t-quark mass or the Z,Z' mixing angle 9m- Although, a direct 
comparison of the form factors of the different approaches may be meaningless. 


The present article contains a precise formulation of the common treatment of standard weak 
loops and of Z' effects together with their realisation in the code ZEFIT0ZFITTER. Since QED 
corrections are model-independent in the sense that they are well-dehned if vector and axial 
vector couplings, mass and width of the Z' are known, the flexible, semi-analytic multi-purpose 
code ZFITTER |]^ can be used after some modihcations for the calculation of QED corrections. 
Thus, we do not discuss explicit formulae on QED corrections in presence of a Z'. 

In section 3, we comment on the Z width. It is shown how the weak form factors which re- 
normalise the Fermi constant {pz) and the weak mixing angle {kz) are influenced by a mixing 
of Z, Z'. In section 4, the same is done for the corresponding form factors of the differential 
cross sections p(s, Q^) and Ke{s,Q‘^), Kf{s,Q^), K,ef{,s,Q‘^). Section 5 contains some explicit 
cross section formulae. In section 6, the structure of the package ZEFIT and its interplay with 
ZFITTER is described. It contains the above mentioned changes of form factors together with 
Born cross section parts containing Z' exchange [for applications at energies beyond the LEPl 
region]. Further, in appendix A the explicit vector and axial vector couplings of the Z' in two 
important classes of extended gauge models, the Eg based and left-right symmetric models, are 
determined as functions of one free parameter, 9e6 or aiR, respectively. Appendix B contains 
the output of a test program for the use of ZEFIT at LEPl energies. 

At the end of this introduction to the subject we would like to stress that the code originally 
is intended for LEPl physics, but may also be used at higher energy, e.g. for Z' searches at 
LEP 200 or LINAC 500. An application of this kind and a comprehensive list of references 


may be found in [11 
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2 Gauge Boson Mixing 

The coupling constants are defined for the symmetry eigenstates Z, Z'. The Lagrangian 




contains currents of the form 


Jn=Yl ['^/(’^) + 'y 5 af{n)] /, n = 7 , Z, Z'. 
/ 


( 4 ) 


(5) 


In the following, complications will arise from a mixing of Z^Z' since our renormalisation 
is performed on mass shell, i.e. for mass eigenstates Z\^Z 2 - In this chapter, we distinguish 
between these states: 

^ Zi \ _ ( cos 6 m sin 6 m \ ( Z 
Z<) I I -sin^M cos 6 m } i Z' 


( 6 ) 


^2 / \ — 

The weak mixing angle 6iv is related to the gauge boson masses and to the gauge boson mixing 
angle 6m as follows: 

M| - Ml 


tan 6m 2 _ j ^2 > 


Mz = 


M, 


w 


These equations correspond to 


Pmix — 




cos 6w 


Ml l + tj^ Ml/Ml 


Ml 


Ml cos^ 6w Ml 1 + ^ V 


(7) 

( 8 ) 

(9) 


The Z boson mass measured at LEPl is Mi = 91.177 GeV. The couplings of the mass 
eigenstates to fermions are : 


g' 

n/(l) = cos 6 'mM + — sin 6*M'y/, 

Vf{2) = cos6MVf - ^sin^MM’ 

g = {V2Gf,Ml)^^‘^, Vf = a/(l - 4|Q/| sin^ 6w), 

with analogue definitions for the axial couplings. 

The photon couplings are defined such that Qe = —1. 


«/ = Pif), 


( 10 ) 

( 11 ) 

( 12 ) 


3 Partial Z widths in Presence of Z, Z’ Mixing 

Without Z, Z' mixing, the matrix element for the decay of the Z boson into a fermion pair 
may be written as follows [§: 




^Mle'^^fttfU 7„(1 + 75)-4sin^6'i4/K/ 


u 




V2 


ze°‘afU 


Vf 

7a75 + 7a — 
«/. 


M, 


(13) 
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where p and k contain the weak loop corrections as determined in the on mass shell renormal¬ 
isation scheme, after the replacement of the coupling constant a by the muon decay constant. 
They are related by: 


2 sin^ 6w cos^ 6w 

The factor (1 — Ar) becomes part of the form factor p [§]. 
The resulting decay width is: 


f/ 


V2l27r 
a/2 bvr 


c/P/ 


c/ 


1 — 4:\Qf \ sin^ 6wi^f + 8(|Q/| sin^ dw^f) 
+ 0^1 , 


2 


where c/ is a color factor in case of quarks. Further, we used that 


^ = 1 - 4|Q/| sin^ Ow^f- 

Qf 

The effective vector and axial vector couplings are dehned as follows: 


h/ = v^/3^(/), 


(14) 


(15) 


(16) 


(17) 




a/ 


1 — 4:\Qf\ sin^ OwKf , 


(18) 


where I^{f) is the weak isospin of fermion /. It has to be stressed here that the weak form 
factors are introduced originally for the overall normalisation of the matrix element and for 
the weak mixing angle, but not for the axial and vector couplings. 


If now Z and Z' mix, the above formulae must be written down for the mass eigenstate Zi 
with mass Mi = 91.177 GeV and in terms of the couplings of the ’physical Z boson’ Zi. One 
can rewrite the Born couplings after mixing as follows: 


®/(^) ~ Cm cif + sm— o,r — [Cm + sm - )o,f, 

9 ^ gaf ' 

a/(l) = (1-!//)«/> 


9 df 9 df 

Uf = —Sm - 1 " (1 ~ Cm) ~ —Sm -• 

gaf gaf 

Similarly, we make the ansatz: 


m(i) 

«/(l) 


Vf + tMVfg'/g 
ttf + tMa'j-g'/g 


1 -4|Q/|sin^ 9w{^ - ^f), 


from which we derive: 




(l-Vf/a,) ‘ 


f vf + tMv'fa'/a 

+ (mo'/JVs 



(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 
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( 24 ) 



Vf{l)=af 1 — AlQflsin^ OwKf . (29) 

The width of the ’physical Z’ in presence of standard weak loops and mixing now is: 

r' r n 

f(l)/ = ^j^Cfpf[l-4\Qf\sm^9^Kf + 8(\Qf\sin^0^Kfr_ 

n /\/f3 

= ^-^c,[h,(l)^ + a,(l)2]. (30) 

This is the expression for a partial width of the Z boson studied at LEPl as they are used in 
the ZEFIT package for hts including the Z, Z' mixing. 

One now can study the expressions which result in case of both small mixing and weak loop 
corrections: 

%(1) = af^ = (1 - yf)^/Ff af = ^fFAcM af + sm «/)• (31) 

A naive ansatz for the same effective coupling could be: 

af^ = CM{^/Pf a/) + SMa'f. (32) 

A direct comparison shows that one should not expect a more than tiny difference. 

For the resulting effective vector coupling, one gets similarly: 

Vf{l) = pf{cMaf + smOj))! - 4sin^ 0^(1 - Xf)K,f]. (33) 

A naive ansatz could be here: 

= CM[y/Wf «/(l - 4sin^6'^y|(5/|fi;/)] + SMv'f. (34) 
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4 Weak Form Factors of the Scattering Matrix Element 
in Presence of Z, Z’ Mixing 


The Born matrix element for the scattering through the mass eigenstate Z\ (being observed 
at LEPl) is: 


Af i(s, cos'd) 


m- 


rae(l)a/(l) 




7/. 


Ae(l) 


+ 75 


7 


i Vfi^) 
V«/(i) 


+ 75 


, (35) 


where ml denotes the complex mass parameter including hnite width effects^. The following 
short notation is used: 

B'^ = [u^A^u^ ■ [ufB'^Uf]. (36) 

A similar ansatz may be written in the t-channel. The matrix element may be rewritten in 
terms of standard theory (unmixed) variables af,j,Av?6w'- 


Af i(s, cos'd) rsj 


, ,2 

-4|Q/| sin^ dw{^ — X 


(1 - ye){l - yf)[L^, ®L^- 4|(5e| sin^diy(l 
f)L^ ® 7^ + lQ\QeQf\ sin'^diy(l - Xe){l - X 




= 7^(1 + 75), 

where again we used that one can write 


^/(l) 

0 /( 1 ) 

0/(1) 


1 -4|Q/|sin^ dw(l 
(l-d/)a/, 


Xf), 


(37) 

(38) 

(39) 

(40) 


with the same dehnitions of x/, ?// as in the case of the Z width. 


Without mixing, the weak loop corrections influence the matrix element as follows: 

Alz(s, cos'd) rsj 

aeafp{s, cos'd)[L^ ® — 4:\Qe \ sin^ 9wKe{s, cos'd) 7 ^ ® 

—4|Q/| sin^ Owi^fis, cos'd)L^ ® 7 ^ + l6\QeQf \ sin^ 9wt^ef{s, cos'd) 7 ^ ® 7 ^]. (41) 

For massless fermions, the four form factors p^K^-iKf, K^f are the most general ansatz for the 
weak radiative corrections. In Born approximation, p = k = 1. 


s — mt 


1^7 


In case of Z, Z' mixing, the matrix element may be written in a similar form: 
A4i(s, cos'd) r\j 


1 


ae(l)a/(l)p^^(s, cosh*) 


s — ml 

[L^ ®L^ -4|Qe|sin^ 9wnf {s, cos ■9)'y^ 




- 4|Q/| sin dvi/K/ (s, cos 79 )L^ ®+ l6\QeQf\sm 9wK^fis, cos'd) 7 ^ ® 7 ^], (42) 


^ We do not discuss here problems connected with the definition of gauge boson masses in dependence on 
the handling of the energy dependence of the width. 
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where it is: 



PmixiX ?/e)(l yf)Pi 

(43) 


= (i-a;/)^/, 

(44) 

^ef 

= (1 -Xe)(l -X/)Ke/. 

(45) 


The matrix element may be rewritten in terms of effective weak neutral vector and axial vector 
couplings: 


Afi(s, cos'd) 


s — nil 




[aea/7^75 0 7^,75 + ® 7^75 + aeV/lpTib ® 

+Ve/7M®7'"], (46) 


a/(l) = v^p“(s, cos'd) (47) 

h/(l) = 4/(1) 1 — 4|Q/| sin^ dwK^(s, cos'd) , (48) 

he/(l) = ae(l)h/(l) + he(l)a/(l) - ae(l)a/(l) 1 - 16|QeQ/| sin'^ cosd) .(49) 

An equivalent notation is (with and without mixing): 


«/ = (50) 

Vf = y^K/'U/+ (1 - K/)a/, (51) 

Vef = pKefVeVf + {Ve- ^/pKefVe)af + tteiVf - ^/pKefVf) - {1 - Kef)aeaf. (52) 


Alternatively, one could dehne the axial vector couplings to be unchanged by the radiative 
corrections. Then, the Fermi constant absorbs the weak form factor p(s, cos'd) and becomes 
dependent on the process and its kinematics: 

^ = G^p^{s, cos^). (53) 

The other form factors renormalise the weak mixing angle sin^ d^v = 1 — 

f sin^ dpi/K^(s, cos'd) 

sin^ Ow I (g,cosd) . (54) 

sin^ 6 w\Jnff{s, cos'd) 


Finally, one should mention that the mixing of gauge bosons Z, Z' as discussed here could be 
mimicked by a mixing of standard fermions with exotic fermions , jl^pl with quite similar 
influence on 


In order to simplify a comparison of the present approach to that of other groups, we now give 
the explicit leading top quark mass m* dependence of the form factors and Ar |T^ |T^ ^ 
in case of mixing: 


sin^ 6wM^ = 


sin^ dw = 


ira 


/{V2G, 


1 — Ar 


Ml 


= 1 


Ml 


Mlpr. 


(55) 

(56) 


^In the common influence of exotic fermions and a Z' is discussed. 






















and pmix is defined in (H). Fnrther|, 


Ar 


1 + 


(1 + 


cos^ 6w 
sin^ 9w 


5p){l - Aa), 


(57) 


6p = 3T 


1 - {271^ 


19)T 


tts 2 _ 
-7 t(1 + 
71 3 



rrij 

a/2 Stt^' 


(58) 

(59) 


Here, as is the strong interaction conpling constant, and Ar'’®"* contains the 0(a) corrections 
to Ar withont the contribntion T. The Aa contains the fermionic one-loop insertions to a. 
For the cross section form factors, analogne formulae hold: 


p^ = 

1 + ApT™ 

Pmixi^ 2/e)(l Vf) ^ y 

(60) 

K^ - 
Kf - 

sm U\y 

(61) 

K^ - 

(1 - x,)(l - X,){1 + A<P)(1 + 

sm 

(62) 

For the partial Z widths. 



Pf 

— Pmix(l Vf) ^ 5 

(63) 


sm Ow 

(64) 


The cross section form factors and those of the partial Z widths differ by their remnant parts 
and, in case of the p, by their mixing factors. At very high energy, the remnant parts of form 
factors become substantial and both the approximate flavor independance, the similarity of 
width and cross section form factors, and the factorization property of Kef are lost. 


We remind the reader here that it was not our intention to discuss the various weak mixing 
angle definitions which allow an aproximate short hand notation for partial Z widths and cross 
sections. We introduced the weak corrections such that they are exact to weak one loop order 
after introduction of the form factors. 

A good approximation to the so-called effective weak mixing angle which often is used for a 
description of LEPl data is 

sm‘^6w,eff = K, (65) 

where one can take any one of the (real part of) form factors k/, calculated at s = M|, cos'd = 0 
or the corresponding form factor from a partial width, e.g. Fg. For some details see and 
for a comparative discussion [^. 


At the end of this section, we should mention that all the above derivations are valid for 
Bhabha scattering 0 ), too. This is due to the fact that the discussion has been based on the 
matrix elements M. before building cross sections out of them. 

^We follow the conventions of ZFITTER4 with the following setting of flags: IAMT4=2, IQCD ^ 0. The 
actual version of ZEFIT is to be used together with ZFITTER v.3.05 until the official release of version 4. The 
corresponding flag setting in version 3 is: IAMT4=1. The resummation of the higher order QCD terms is not 
realized there. 


9 












5 Improved Born Approximation in Presence of Gauge 
Boson Mixing 

The matrix element M.z' with Z' exchange is in case of mixing: 

/2 

M2(s,cosi9) ~ - 2 [®e(2)a/(2)7^75 0 7^75+ he(2)a/(2)7^ 0 7^75 

5 — ^^2 

+ ae(2)h/(2)7^75 07 '" + Vef( 2 )j^ ® 7 ^"], ( 66 ) 

where a/(2),h/(2) are the renormalized vector and axial vector couplings of the Z 2 . 

As long as we search for a Z' well below the production threshold, we can safely neglect here 
all radiative corrections to this amplitude, 

o;(2)=o;(2), ty(2)^v,(2), v,f(2) = v,(2)v,(2). (67) 


In sum, our discussion leads to the following net matrix element, where we also add up the 
photon exchange diagram with running QED coupling: 

M = M^ + Mi + M2. (68) 


Let us now have a look at the squared matrix elements. In case of massless fermion produc¬ 
tion, four different combinations of coupling constants may occur in reactions (|I],^) from the 
interference of the vector bosons m and n: 


Ci{m,n) = ae{m)al{n)af{m)a*jr{n) + ae{m)al{n)vf{m)v*jr{n) ( 69 ) 

Ve{m)v*{n)af{m)a*j{n) + Vef{m)vlf{n), 

C2{m,n) = ae{m)v*{n)vf{m)d*f{n) + Ve{'nT-)al{n)af{m)Vf{n) ( 70 ) 

-f de{m)vlf{n)df{m) + a*(n)he/(m)a)-(n), 

Cz{m,n) = de{m)v*{n)df{m)d*jr{n)+Ve{m)dl{n)df{m)d*f{n) ( 71 ) 

de{m)v*f{n)vf{m) + dl{n)vef{m)vj{n), 

C/^{m,n) = de{m)dl{n)df{m)vf{n) + de{m)dl{n)vf{m)d*j{n) ( 72 ) 

- 1 - Ve{m)v*f{n)df{m) + vl{n)vef{m)d*f{n). 


The starred couplings of vector boson n and the corresponding propagators are complex con¬ 
jugated; a procedure which is necessary only in the s-channel. 

Now we have all formulae needed to write down the improved Born cross section. For 
initial state radiation, they read^: 


a 


0 

T 


TTCK 


N 


2s 


-3fJe Xm{s')Xn{s') Ci{m,n)LiHi 


m,n=0 


+ Cs{m, n)L2Hi + C^{rn, n)LiH2 + C2{m, n)L2H2 


( 73 ) 


While for final state radiation the s' is replaced by s, the s-dependence is more complicated for the 
initial-final state interference We further remind the reader that ZFITTER doesn’t return differential 

cross sections at all; (^ is shown for illustrational purposes. 
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(74) 


a 


0 

FB 


Xn('S) 

da 

d cos’d 


Fa 


N 


2s 


-3fJe Xm{s')Xn{s') C 2 {m,n)LiHi 


m,n=0 


+ Ci{m, n)L2Hi + C3(m, n)LiH2 + C'i(m, n)L2H2 , 

al s 

4Fa s — rn^ ’ 


(1 + cos^'i9)(Ty + (2 cos'(9) cr%. 


(75) 

(76) 


The quantities Li, L 2 , Hi and H 2 are combinations of the polarizations of the beams (e , e"*") 
and of the helicities of the hnal fermions: 


Li — \ — A_|_A_, L 2 — A_|_ — A_, 


= j(i 


h+h.), 


H, = -(Ah 


A-), 


(77) 


with A_ and A+ {h_ and h+) being the polarizations of electron and positron (fermion and 
antifermion). The masses in the propagator Xn('S) are: 

mg = 0 7 

ml = Mf - iMiTi Z 

ml = M| — iM2T2 Z'. 


In case of a {Z,Z) Born cross section with a polarization of the electron beam, the following 
combinations correspond to the well-known coupling factors: 

Ci{Z, Z) + X+Cs{Z, Z) = {al + vl + 2\+a,v,){a) + nj), (78) 

C2{Z, Z) + Xj^Ci{Z, Z) = pOeUe -|- A_|_(ag -1- vl)]{2afVf). (79) 

With the above dehnitions, we have all the necessary prerogatives to calculate cross sec¬ 
tions with both weak loop effects and Z' exchange. The resulting Born formulae can be used 
as input for a QED calculation. 


The above expressions for and a^^ are used in the redehnitions to be performed in the 
subroutine BORN of ZFITTER. 


6 Structure of ZEFIT 

The package ZEFIT should be run together with ZFITTER. For execution it has to be loaded 
before the package ZFITTER since it contains subroutines which originate from ZFITTER but 
are modihed for the description of the Z'. 

Searching for a Z', we assume that the Z boson interactions are correctly described by the 
Standard Model. Then there are two ways to use the package. Either one performs a model- 
independent ht to data with ZFITTER and searches with the result - partial widths or effective 
couplings - for a Z', using basically the ZWRATE or ROKANC subroutines of DIZElf], or one 
tries an immediate ht to the cross sections, using basically subroutine ZCUT calling EWINIT 

® ZWRATE calls ROKAPP to calculate the weak form factors of the widths. 
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of ZFITTER|^. The common package ZEFIT0ZFITTER is prepared for both applications. The 
second approach, within the (possibly extended) Standard Model terminology, is technically 
more involved. So, we will use for the present purpose of demonstration only the correspond¬ 
ing interface and branch with subroutine ZANALY of ZFITTER v.3.05. 

Input arguments: 


AMZ is the mass of the Z boson, 

AMH is the mass of the Higgs boson in GeV, 

AMT is the top mass in GeV, 

INDF selects the hnal state fermion type, see also the table in the test example in Ap¬ 
pendix B, 

SQS is the centre of mass energy in GeV, 

IFAST = 1 allows a fast calculation without geometrical and kinematical cuts. It must be 
IFAST = 0 if cuts are required, 

IRCUT = 0,1 chooses between an acollinearity cut and a cut on the photon energy. 

New additional input parameters coming from a Z'\ 

AMZE is the mass of the Z' in GeV, 

ZMIX is the mixing angle between Z and Z', 

IZE = 0,1,2 chooses the model (Standard Model, Eg model, left-right model). 

As examples, we have foreseen a Z' coming from an Eg-GUT and from a left-right symmetric 
model. The model chosen must be specihed by the flag IZE. The Standard Model is realised 
for IZE = 0, the Eg model is chosen with IZE = 1 and the LR model with IZE = 2. Inside 
the Eg model the parameter TETAE6 must be set as the mixing angle of the two extra Z 
generators. For the left-right symmetric model the parameter ANGLR has to be dehned. 
ANGLR is limited, i.e. < ANGLR < a/2. The unit for the parameters TETAE6 and 
ANGLR in ZEFIT is radian. 


Additioual flags of ZFITTER The flags for weak loop and QED corrections are 
subroutine ZINITF 0: 


set m tne 


IWEAK = 0 or 1, switches the 0{a) weak loops. 

IHVP = 1, 2, 3 characterises the vacuum polarization parametrization. The best choice is 

IHVP = 3. 

IQCD = 0, 1, 2, 3, 4 gives the QGD corrections to the vector boson self energies in Ar, the 
widths and the cross section. IQCD = 3 is recommended for LEP applications. 

® EWINIT uses the subroutine ROKAP in ZFITTER v.3.05 (RDKANC in ZFITTER v.4) to calculate the form factors 
of the cross sections. 
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IAMT4 = 0, 1; the leading two loop effects of the type can be included. 

IBOX = 0, 1; WW and ZZ box corrections may be taken into account. 

IFINAL = 0, 1 chooses between the approximate hnal state correction by the factor [1 + 
37r(5//(47r)] and the correct 0{a) hnal state correction including soft photon exponenti¬ 
ation. 

INTERF = 0, 1 switches the 0{a) initial - hnal state interference. 

IPHOT2 = 0, 1; second order QED leading logs may be taken into account. 

Additionally to these hags of ZFITTER, we introduced 

lALFRUN = 0,1; which allows an independent switching oh and on of running aqEo- 


Output parameters: 


SBORN - Born cross section in nb. 

STOT - total cross section in nb. 

ABORN - Born forward-backward asymmetry. 

ATOT - forward-backward asymmetry. 

In hg. 1 the block diagram of ZEFIT is shown. The subroutines in the dashed boxes remained 
unchanged. 
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|ewcoup| 


I SETFUNI 

! SFAST ! ! SCUT ! 
BORN 


Figure 1: Structur of the package ZEFIT. Subroutines in solid boxes are rewritten and 
replace the corresponding subroutines of ZFITTER in case of Z, Z' mixing. The subroutines in 
dashed boxes remain unchanged. They are shown here in order to illustrate the interplay of 
ZEFIT with ZFITTER. 
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Appendix A 

Additional parameters in theories with extra Z bosons 

The couplings of Z\ and to fermions are in general a result of a mixing between the Standard 
Model Z and the Z' as introduced in (|): 

Vf{l) = cosOm Vf + sin6 'm —v'f, 

g' 

0/(1) = cosOm + sin^jvr —«/, 

Vf{2) = cos6m v'f - siuOm ^Vf, 


0 /( 2 ) = cos 6 m o!j — sindM 


(A.l) 


The couplings of the extra Z boson to fermions O/ and v'j depend on the particular model. 
In ZEFIT we have implemented extra Z bosons coming from an Eq GUT or from a left-right 
model. 


In the Eq gut the Eq group M] is assumed to be broken to the standard model group 


structure in the following way |]W 


Ue —^ ^O(IO) X U{1)^ —^ SU{5) x U{1)^ x U{1)^ —^ 
— ^ SU{3)c X SU{2)l X U{1)y x U{1)^ x f/(l)^. 


(A.2) 


We further assume that the following linear combination of the two extra Z bosons Z^ of 
U{1)^ and Z^ of U{1)^ is light and can be detected at future colliders [p^l : 


Z' = cosOe Z^ -|- sin^E Z, 


E ^ip- 


(A.3) 


The couplings between the Z' and fermion / are K 


< = 

^Q^cosOe + |Q-i,sin6', 

K = 

|Q^COs6'£; -F 

< = 

Q^cosOe + Q^rsin^, 

< = 

2Qy^ cos 9e, 

< = 

—Qj^cosOe + Q^rsin^, 

K = 

0, 

a'd = 

Q^cosOe + Q^sin^, 

v'd = 

—2Q^ cos 9e, 


1 


1 



Qxjj 

Te' 


(A.4) 


(A.5) 

The Z' coupling constant g' may be determined by the assumption that the renormalisation 
group evolution of g' is the same as that of g p2| , 


^' = y 2 sin Ow g- 


(A.6) 
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Of particular interest for applications is the completely specified case Z' = — 

Y^5/8 Z^ as suggested by superstring theories |^, The most general Z' in a Eq GUT 
can therefore be described by three additional parameters: Om =ZMIX, 6e =TETAE6 and 
M(Z2) =AMZE. 

Another origin of a Z' could be a left-right symmetric model 


25[|. As in the Eg case, we have 


6m =ZMIX and M(Z 2 ) =AMZE and one additional parameter a =ANGLR. In this model we 
have the following couplings of the Z' to the fermions: 


< = 

1 





K = 



a'e = 

1 

2 “’ 

< = 

1 

a 


< = 

1 

—a, 

2 ’ 

< = 

1 

3a 

- 


1 


1 

/ 

a'd = 

2 “’ 

v'd = 

3a 



The relation between g and g' is: 

g' = sin Ow g. 


(A.7) 


(A.8) 


Goncerning the width of the Z', we will assume that it can decay only into particles of the 
known three fermion generations including the top-quark. 


Of course, one also can make no model assumptions about the Z' at all. Then, one has to 
specify arbitrarily all needed couplings of the Z' to fermions and the width of the Z'. In 
addition, the Z, Z' mixing angle and the Z' mass have to be given. 
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Appendix B 

The Main Routine of the Test Program 

PROGRAM ZEFIT 
c 
c 


** 



** 

** 

ZEFIT 

a package for extra Z searches at LEP 

** 

** 



** 

** 



** 

** 

A package 

[1] ZPRIME physics in presence of QED and weak 

** 

** 

loop corrections as realised in ZFITTER [2] 

** 

** 

A possible 

mixing of Z and Z’ is taken into account. 

** 

** 



** 

** 

Version: 

2.00 

** 

** 

Release: 

November 25,1991 

** 

** 

Authors: 

A. Leike, S. Riemann, T. Riemann 

** 

** 


Inst, for High Energy Physics, Zeuthen 

** 

** 

Contact in 

L3: S. Riemann, riemannsOcernvm 

** 

** 



** 

** 

[1] A.Leike, S.Riemann, T. Riemann: 

** 

** 

Z.Z' 

Mixing 

** 

** 

in Presence of Standard Weak Loop Corrections, 

** 

** 

Munich Univ. prepr. LMU-91/06 (Dec 1991) 

** 

** 

[2] D.Bardin et.al.: 

** 

** 

A Users Guide to ZFITTER: An Analytical Program for 

** 

** 

Fermion Pair Production in e~+ e"- Annihilation, 

** 

** 

in preparation. 

** 

** 



** 


* 

IMPLICIT REAL*8(A-H,0-W,Z) 

IMPLICIT COMPLEX * 16 (X,Y) 

REAL*4 TIME1,TIME2 
PARAMETER (NSARR=10,NM0DEL=1) 

* 

* NSARR - NUMBER OF SQRTS POINTS TO BE CALCULATED 

* NMODEL - NUMBER OF POINTS SPECIFIED WITHIN AN EXTRA Z MODEL 

* 

C0MM0N/R0VEFZ/ARR0FZ(0:10),ARVEFZ(0:10) 

COMMON/FRINIT/ NPAR(30),ZPAR(30) 

COMMON/ZE/TETAMD,AMZE,GAMZP,ZMIX,CHIE 
+, VEEZE,VEFZE,AEFZE,VQEZE(6),VQFZE(6) 

+, XVEEZE,XVEFZE,XAEFZE,XVEFGZ,XAEFGZ 

+, XVQEZE(6),XVQFZE(6),XVQFGZ(6) 
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COMMON /IZE/ IZE 

* ARR0FZ(0:9) - CALCULATED EFFECTIVE RO'S FOR EACH CHANNEL 

* ARVEVZ(0:9) - CALCULATED EFFECTIVE VECTOR COUPLINGS FOR EACH CHANNEL 

* ARR0FZ(0:10) AND ARVEFZ(0:10) ARE UNDEFINED AND NOT USED 

* ARVEVZ = 1 - 4 * ABS(QF) * SINTW**2 ( SINTW**2 IS EFFECTIVE ) 

DIMENSION WIDTHS(0:10) 

DIMENSION TETAE6(NM0DEL),ANGLR(NMODEL) 

DIMENSION SARR(NSARR) 

* ARRAY OF CMS ENERGIES 

DATA SARR /87.dO,88.DO,89.ODO,90.dO,91.dO,91.17d0,92.dO,93.dO 
+, 94.d0,95.d0/ 

** 

* CHOICE OF E_6 ZPRIME MODEL (TETA6 in RADIAN) 

C DATA TETAE6 /O.DO, 1.5708D0, -.91174D0/ 

DATA TETAE6 /O.DO/ 

* CHOICE OF LEFT-RIGHT EXTRA Z MODEL (ANGLR in RADIAN) 

DATA ANGLR /O.81649658D0, l.DO, 1.4142135D0/ 

C sqrt(2/3) sqrt(2) 

** 

* INPUT FOR INITIALIZATION SUBROUTINE ZINITF: 

* Z-MASS (ALL INPUT MASSES IN GEV) 

AMZ=91.180D0 

AMZLEP=AMZ 

* TOP-MASS 

AMT=150D0 

* HIGGS-MASS 

AMH=3D2 

* ZPRIME MASS (NOT RELEVANT IF IZE = 0) 

AMZE= 5000.DO 

* TWO QCD-CORRECTION FACTORS (THE SECOND ONE FOR B-BBAR CHANNEL) 

QCDC0R=1.040D0 

QCDC0B=1.045D0 

* 

* IZE = 1,2 (0) EXTRA Z CONTRIBUTIONS ARE (NOT) INCLUDED 

* =0 STANDARD MODEL 

* = 1 E_6 MODEL (TETAE6 MUST BE SPECIFIED) 

* =2 LEFT-RIGHT MODEL (ANGLR MUST BE SPECIFIED) 

IZE=1 

* ZMIX - Z, ZPRIME MIXING ANGLE 

ZMIX = O.OIDO 

* IF IFAST=1 WITHOUT CUTS 

* IF IFAST=0 CUTS ARE POSSIBLE 

IFAST=1 

* PLEASE, LOOK INTO ZINITF IN ORDER TO INITIALIZE ALL FLAGS AND 

* PARAMETERS WHICH COULD BE (AND SHOULD BE) INITIALIZED FOR A SPECIFIC 

* TASK. 

* 
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IF(IZE.EQ.O) THEN 
INMOD=l 
ELSE 

INMOD=NMODEL 
END I 

* SPECIFY THE MODEL DEPENDENT ANGLE 

DO 20 IZP=1,INM0D 

* EXTRA Z FROM E_6 MODEL 

IF(IZE.EQ.l) TETAMD=TETAE6(IZP) 

* EXTRA Z FROM LEFT-RIGHT MODEL 

IF(IZE.EQ.2) THEN 

ANGMIN = DSQRT(2.D0/3.D0)-l.D-5 
ANGMAX = DSQRT(2.D0)+l.D-5 
TETAMD=ANGLR(IZP) 

IF(TETAMD.lt.ANGMIN .OR. TETAMD.GT.ANGMAX) THEN 
PRINT*,’ LEFT-RIGHT PARAMETER OUT OF RANGE’ 

STOP 
END IF 
END IF 

* 

CALL ZINITF(AMZ,AMT,AMH,QCDCOR,QCDCOB,SW2,WIDTHS) 

* OUTPUT OF ZINITF: 

* SW2 - CALCULATED QUANTITY = 1-AMW**2/AMZ**2 

* WIDTHS(INDF) - CALCULATED PARTIAL CHANNEL Z-WITDHS (IN MEV) 

* FOR 0<INDF<9 AS DESCRIBED ABOVE 

* FOR INDF=10 IT CONTAINS CALCULATED TOTAL Z-WIDTH 

* 

IF(IZE.EQ.O)PRINT’(/, ” STANDARD MODEL, NO EXTRA Z EXTENSION”,/)’ 
IFdZE.EQ. 1) PRINT’(/, ” EXTRA Z FROM E6 GUT”,/) 

IF(IZE.EQ.2)PRINT’(/, ” EXTRA Z FROM LEFT-RIGHT MODEL”,/)’ 

PRINT 1000,AMZLEP,AMT,AMH 

IFdZE.EQ. 1) PRINT 1100,AMZE,ZMIX,TETAMD 

IFdZE.EQ. 2) PRINT 1200,AMZE,ZMIX,TETAMD 

PRINT 1300,QCDCOR,QCDCOB,SW2 

PRINT 1500,WIDTHS 

PRINT 1550,NPAR(1),NPAR(2),NPAR( 3),NPAR( 4) 

&, NPAR(8),NPAR(9),NPAR(10),NPAR(16) 

* 

* GEOMETRICAL CUTS OVER THE ANGLE BETWEEN E+ AND OUTGOING ANTI-FERMION 

ANG1=140D0 

ANG2=40D0 

IF (IFAST .NE. 1) PRINT 1700,ANG2,ANG1 

* 

* KINEMATICAL CUTS: 

* IF IRCUT=1 THEN SPRIME (IN GEV**2) 

* IF IRCUT=0 THEN ACOL+EMIN (IN DEGREES AND GEV, RESPECTIVELY) 

* 
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IRCUT = 0 
SPRIME= lD-6 
ACOL = 25.DO 
EMIN = .ODO 
IF(IFAST.EQ.O) THEN 
IF(IRCUT.EQ.l) THEN 
PRINT 1800,SPRIME 
ELSE 

PRINT 1900,ACOL,EMIN 
END IF 

ELSEIF(IFAST.EQ.l) THEN 
PRINT 1600 
ELSE 

IFAST = 1 
PRINT 1600 
END IF 

* 

DO 11 INDF=2,2 


* 

INDF=0 

FOR 

NEUTRINO 

* 

=1 

FOR 

ELECTRON 

* 

=2 

FOR 

MUON 

* 

=S 

FOR 

TAU 

* 

=4 

FOR 

UP 

* 

=5 

FOR 

DOWN 

* 

=6 

FOR 

CHARM 

* 

=7 

FOR 

STRANGE 

* 

=8 

FOR 

TOP ( RETURNS ZERO ) 

* 

=9 

FOR 

BOTTOM 

* 

=10 

FOR 

HADRONS 


PRINT 2000,INDF 
GAMEE=WIDTHS(1)/IDS 
GAMZ=WIDTHS(10)/1D3 
IF(INDF.NE.10) THEN 

GAMFI=WIDTHS(INDF)/IDS 
ELSE 

GAMFI=ODO 
DO 4 IH=4,9 

4 GAMFI=GAMFI+WIDTHS(IH)/1DS 

END IF 

IF(INDF .NE. 10) PRINT 2500 
IF(INDF .EQ. 10) PRINT 2600 

* DO LOOP OVER THE SQRT(S)=SQS 

DO 10 IS=1,NSARR 
SQS=SARR(IS) 

* 

* NOW EVERYTHING WILL BE CALCULATED IN THE FRAMEWORK OF THE STANDARD 

* MODEL BY ZANALY- SUBROUTINE. 


22 



IF(INDF.NE.IO) THEN 

CALL ZANALY(IFAST,INDF,SQS,IRCUT,SPRIME,ACOL,EMIN,ANGl,ANG2, 

& AMZ,AMT,AMH,QCDCOR,QCDCOB,SHORN,STOT,ABORN,ATOT) 

CS =STOT 
AFB=ATOT 

IF(INDF.EQ.O) THEN 
PRINT 3000,SQS,CS 
ELSE 

PRINT 3100,SQS,CS,AFB 
END IF 
ELSE 

CS =0D0 

DQ 53 INDFH=4,9 

CALL ZANALY(IFAST,INDFH,SQS,IRCUT,SPRIME,ACOL,EMIN,ANGl,ANG2, 
& AMZ,AMT,AMH,QCDCOR,QCDCOB,SHORN,STOT,ABORN,ATOT) 

CS =CS +STOT 
53 CONTINUE 

PRINT 3000,SQS,CS 
END IF 

10 CONTINUE 

11 CONTINUE 
20 CONTINUE 

1000 FORMATClX,/IX,’INPUT : MZ = ’,F6.3,’ GEV’ 


& 


MT 

II 

2,’ GEV’ 





MH = 

’,F8.2, 

’ GEV’ 



1100 FORMAT 

(’ 


AMZE 

= ’,F6.1,’ 

GeV’ 


&, 


) 

ZMIX 

= ’,F6.4,’ 

TETAE6 = ’ 

,F8.4,/) 

1200 FORMAT 

(’ 


AMZE 

= ’,F6.1,’ 

GeV’ 


&, 


) 

ZMIX 

= ’,F6.4,’ 

ANGLR = ’ 

,F8.4,/) 

1300 FORMAT 

(’ 

QCD CORRECTION 

FACTORS = 

’ ,2(F7.5,2X),/ 


& 

> 

OUTPUT: 

SW2 

= 

’,F5.4) 



1500 FORMAT (IX,’PARTIAL AND TOTAL Z-WIDTHS IN MEV’,/ 

&,5x,’nu,nubar =’,F7.l,8x,’e+,e- =’,F7.1,6x,’mu,mubar =’,F7.1 

&,/,5x,’tau+,tau- =’,F7.l,8x,’u,ubar =’,F7.l,8x,’d,dbar =’,F7.1 
&,/,5x,’c,cbar =’,F7.l,8x,’s,sbar =’,F7.1 

&,/,5x,’t,tbar =’,F7.l,8x,’b,bbar =’,F7.l,8x,’TOTAL =’,F7.1 

&,/) 

1550 FORMAT (/,’ FLAGS IN INITIALISATION ROUTINE ZINITF:’, 

& /,’ IWEAK =’,I2,’ IHVP =’,I2,’ IQCD =’,I2, 

& ’ IBOX =’,I2,/,’ INTERF =’,12,’ IFINAL=’,I2, 

& ’ IPH0T2 =’,12,’ IAMT4 =’,I2) 

1700 FORMAT (IX,’GEOMETRICAL ACCEPTANCE CUTS: ANGl,ANG2=’,2(F5.1,2X)) 
2500 FORMAT (/,IX,’SQRT(S) [GeV]’,6X,’CROSS SECTION [nb]’,4X, 

& ’ASYMMETRY’) 

2600 FORMAT (/,IX,’SQRT(S) [GeV]’,6X,’CROSS SECTION [nb]’ ) 

1800 FORMAT (/,IX,’KINEMATICAL CUT: SPRIME=’,E6.1,IX,’GEV’) 
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1900 FORMAT (/,IX,^KINEMATICAL CUTS: ACOL=^F5.1,IX,’DEGREES2X, 
& ’EMIN=’,F5.1,1X,’GEV’) 

1600 FORMAT (/,1X,’N0 KINEMATICAL CUT (IFAST=1)’) 

2000 FORMAT (IX,/,’ CHOSEN FERMION CHANNEL INDEX - INDF=’,I2) 

3000 FORMAT (7X,F7.3,IIX, F13.4) 

3100 FORMAT (7X,F7.3,11X,2F13.4) 

END 


Output of the Test Example 


EXTRA Z FROM E6 GUT 


INPUT : MZ = 91.180 GEV MT = 150.00 GEV MH = 300.00 GEV 

AMZE = 500.0 GEV ZMIX = 0.0100 TETAE6 = 0.0000 

QCD CORRECTION FACTORS = 1.04000 1.04500 

OUTPUT: SW2 = .2260 

PARTIAL AND TOTAL Z-WIDTHS IN MEV 


nu,nubar 

169.2 

(D 

+ 

CD 

1 

II 

83.2 

mu,mubar = 

83.2 

tau+,tau- = 

83.0 

u,ubar = 

296.5 

d,dbar = 

386.0 

c,cbar 

296.1 

s,sbar = 

386.0 



t,tbar 

0.0 

b,bbar = 

379.8 

TOTAL = 

2501.6 


FLAGS IN INITIALISATION ROUTINE ZINITF: 

IWEAK = 1 IHVP = 3 IQCD = 3 IBOX = 1 
INTERF = 1 IFINAL = 1 IPH0T2 = 1 IAMT4 = 1 

NO KINEMATICAL CUT (IFAST=1) 


CHOSEN FERMION CHANNEL INDEX - INDF= 2 


SQRT(S) [GeV] 

CROSS SECTION [nb] 

ASYMMETRY 

87.000 

0.1393 

-0.3713 

88.000 

0.2109 

-0.2950 

89.000 

0.3671 

-0.2071 

90.000 

0.7519 

-0.1122 

91.000 

1.4042 

-0.0192 

91.170 

1.4533 

-0.0047 

92.000 

1.1713 

0.0566 

93.000 

0.6813 

0.1090 

94.000 

0.4338 

0.1448 

95.000 

0.3076 

0.1703 
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